Physics of the RVB (Pseudogap) State in the Doped Mott Insulator: Spin-Charge 

Locking 
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The properties of the pseudogap phase above T c of the high-T c cuprate superconductors are 
described by showing that the Anderson-Nambu 57/(2) spinors of an RVB spin gap "lock" to those of 
the electron charge system because of the resulting improvement of kinetic energy. This enormously 
extends the range of the vortex liquid state in these materials. As a result it is not clear that the 
spinons are ever truly deconfmed. A heuristic description of the electrodynamics of this pseudogap- 
vortex liquid state is proposed. 



The hole-doped cuprate superconductors have, over 
the past decade, been shown to exhibit an unusual phase 
in the region of their phase diagram immediately above 
the superconducting "dome" in temperature, and espe- 
cially for underdoped compositions. This region has been 
dubbed variously the "pseudogap" , "spin gap" , and (by 
those unafraid of committing themselves) the "RVB" or 
resonating valence bond region. Fig. 1 shows a gen- 
eralized phase diagram summarizing the fairly generally 
accepted phase boundaries between regions of the com- 
plex behaviors exhibited by these fascinating systems, 
and Fig. 2 (due to Ong and Wangi) is an experimental 
delineation of the regions in which this phase exhibits 
most clearly the behaviors I will discuss later, namely 
a pronounced vortex Nernst effect and nonlinear dia- 
magnetic susceptibility, and momentum-dependent en- 
ergy gap around the Fermi level, in addition to the pseu- 
dogap in density of states and magnetic gap known to 
exist over the wider region below the "pseudogap tem- 
perature" T*. 

In a sense, the existence of a phase of some such sort 
was predicted quite long before it was observed, in that 
it was postulated a few months after the discovery of 
the high Tc cuprates that a quantum liquid of singlet 
pairs of electrons ("valence bonds") might exist in an 
S = 1/2 Mott insulating system, from which the high T c 
superconductor would develop via doping 2 . Baskaran et 
a&, and Fukuyamai even drew schematic phase diagrams 
for such a system which remotely resembled Fig 1. In a 
remarkable and prescient paper—, Kotliar and Liu, using 
slave boson methods, actually predicted the existence of 
a pseudogap phase very much resembling that derived 
by simpler methods in the present work . But of course, 
they could not compare their results with experiment. 

The existence of the pseudogap was first seen in in- 
frared data but not recognized*, and demonstrated by 
magnetic^, tunneling^, thermal^, ARPES 10 and other ex- 
periments over the years. Here I particularly want to 
bring out two observations: first the data of Ong et al on 
the giant Nernst effect and the nonlinear diamagnetic 
susceptibility^! which demonstrate that over a consider- 
able region this phase exhibits a kind of supercurrent in 
a magnetic field; and second the data of Timuslii 2 " which 
show a marked lengthening of the conductivity relaxation 



time. 

The theory I am going to present for this state starts 
from the RMFT (renormalized mean field theory) of the 
superconducting ground state and uses this as an effective 
Ginzburg-Landau energy controlling the fluctuations of 
the variational and other parameters of the state. The 
variational "gap equations" which determine the energies 
of the excitations and the energy of the state areiii4 
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Here, Ek is the quasiparticle energy, its band- 
theoretical kinetic energy referred to the Fermi energy, 
g is the kinetic energy renormalization factor due to the 
modified occupancies (g = 2x/(l + x), x the doping), and 
gj the renormalization factor for J. 

The physics of strongly-correlated electrons is plagued 
by the tendency of limiting cases to exhibit special be- 
haviors. The present case is confused by the fact that the 
antifcrromagnetic insulator is marginally lower in energy 
than the RVB one, at exactly zero for g. Nonetheless we 
choose to ignore that fact and treat the problem using 
the equations Q i n the approximation j«0, since the 
kinetic energy gained in the superconducting-RVB state 
for g finite soon outweighs this small difference. 

The lowest energy translationally invariant state is 
given by these equations with ^ oc cos k x + cos k y and 
Afc cx cos k x — cos ky or vice versa. As was emphasized in 
the early work 2 ^, the solutions, and even the equations, 
are not unique. It is most evident that <; and A may be 
interchanged, or in fact any orthogonal transformation 
on the two is allowed. (It is, however, only in the nearest 
neighbor square lattice without doping that they need 
have the same magnitude.) But this is only a symptom 
of the gigantic symmetry of the solutions, which follows 
from the fact that the Fermion representation on which 
they are based is enormously overcomplete. This over- 
completeness generates a local SU (2) symmetry 2 ., which 
follows from the fact that a spin up hole and a spin down 
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electron create the same state of the insulating Heisen- 
berg model in Gutzwiller projection. This symmetry is 
a true gauge symmetry, in that the states generated by 
it are physically the same state. For translationally in- 
variant, singlet BCS states like JTJ, the gauge symmetry 
implies symmetry under SU (2) rotations operating on 
the three Anderson-Nambu pseudospin vectors t\, t%, T3, 
multiplying respectively the real and imaginary anoma- 
lous self-energies and the kinetic energies of a gas of 
"spinons" . 

But as I said, this is only a tiny fraction of the avail- 
able representations. Also equivalent are states describ- 
able as a staggered flux phase, a true 7r-flux phase, a 
d-density wave, etc., if one allows translationally non- 
invariant solutions or T non-invariant ones. What was 
striking about all good RVB solutions was that they ex- 
hibited four nodes at (tt/2,tt/2) in the spinon energy as 
a function of momentum, and this is simply because of 
their equivalence. The physics behind this fact is that in 
order to describe a projected state each self-energy must 
average to zero hence have a line node; and that the opti- 
mum BCS solution must have as little variation of | A| 2 as 
possible, hence A should be the sum of orthogonal func- 
tions. Without being able to supply a mathematician- 
style proof, this seems to tell us that RVB's will always 
have nodes in their spectra, even though there is no 
Goldstone-type argument for nodes in terms of broken 
symmetry. 

Now we allow g to become finite. Our approach is to 
set up a second triad of r 's, thinking if you like of g 
as defining the T3 direction in this second space. In the 
RMFT ground state wave-function written as^, 

*G=n(^+^ c i T ct -feii o > . ( 2 ) 

k 

the phase <E> of the pair wave-function is normally taken 
to be the phase oiv/u. Our procedure will be to ascribe 
the phase to g and leave u and v real. We can think 
of g as a vector in the t' space which rotates around in 
the t[,t 2 plane, thereby defining the third direction. Its 
meaning at least in the wave function seems to be 
that it is the hole pair amplitude, which in the ground 
state is Bose condensed. 

We now consider the two triplets of vectors r and r' 
to be oriented at arbitrary Eulerian angles with respect 
to each other in the abstract space. The selfconsistent 
solutions for <; and A depend only weakly on these an- 
gles for small g, so a good approximation to the effective 
Hamiltonian for quasiparticles is then 

H = ^ \9 e k + ( T 3 ■ T 2)^k + (73 • n)A fc ] (n k + n-k - 1) 
k 

+ [(n • r[)A k + {t{ • r 2 ) ft ] 4 T cL fci + h.c. 

(3) 

Here we are using a kind of Eilenberger method where 
we derive an effective free energy as a functional of the 



parameters in the one-electron Hamiltonian which leads 
to the gap equation-the most important of these param- 
eters being the angles between the two triads of vectors. 
The "gaps" themselves, ? and A , do not vary much 
at small g, though A, which will become the supercon- 
ducting gap, will of course decrease appreciably toward 
optimal doping. The interaction term in the Eilenberger 
energy, which will be of the form (V 2 + A 2 )/ J , will clearly 
not change much as a function of angles, so that mostly 
we have to consider the effective one-quasiparticle energy. 
The "extended s-wave" solution ^ to the RVB gap equa- 
tions strongly resembles the one-particle kinetic energy- 
which is not terribly surprising since J is a second-order 
consequence of the kinetic energy, and therefore the in- 
teraction term naturally factorizes into terms resembling 
the kinetic energy. (What may be more germane to the 
special nature of the cuprates is the existence of the sec- 
ond, orthogonal, <i-wave solution.) Therefore it becomes 
quite obvious that when t 2 |j 73, the kinetic energy is very 
considerably enhanced, because the coefficient which 
is the multiplier of the first term in (|3J), then becomes 
the sum of two terms of the same structure and the same 
sign. (This was of course the choice which was made 
in the original paper on the RMFT 13 .) Therefore with 
this choice we maximize the quasiparticle energies, which 
necessarily minimizes the total energy. The "locking" is 
illustrated in Fig. 3. 

This is the central result of the theory: that upon dop- 
ing the t 2 axis of the RVB triad locks to the t' z axis of 
the Nambu triad of the charge degrees of freedom of the 
real electrons, so that the T\ axis necessarily lies in the 
t[, t 2 plane and A serves as a ( necessarily real) anoma- 
lous self-energy for the actual electrons, which live on the 
unprimed triad. The reason for the name "spin-charge 
locking" is the similarity to the mechanism of "color- 
flavor locking" of Wilczeki&, which in turn is based on 
the "spin-orbit" locking mechanisms of Leggett for He- 
3i£. But the analogy cannot be carried too far-the pecu- 
liar gauge symmetry of the RVB is unique. 

A second scale is the energy which maintains the bro- 
ken true phase symmetry. The phase in the t[ , t 2 plane 
is meaningful only as a relative variable and acquires no 
stiffness from J- a fact which is equally true for conven- 
tional superconductors, for which the stiffness p s (V0) 2 
is independent of gap parameters or interactions. The 
loss of this stiffness may formally be seen as an "unlock- 
ing" process also, but is not in principle different from a 
conventional T c . There are therefore two transitions and 
three different phases: the low temperature phase is the 
true superconductor where the two triads may be thought 
of as locked together; this undergoes loss of phase coher- 
ence at Tc, but the T3 and t 2 directions remain locked 
together in the spin-charge locked pseudogap phase, so 
that the fluctuating anomalous self-energy remains of d 
symmetry and has nodes along the (n, tt) line; and finally, 
complete unlocking where there may be a pseudogap but 
it does not show momentum dependence-there is no con- 
nection between the spin gap structure and the hole or 
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particle nature of the excitations. It would be only in 
this unlocked phase that we should think of dcconfined 
spinons. 

The transition temperature for the "unlocking" tran- 
sition may be estimated from the energy involved and 
seems unlikely to be sharp. It must lie between T c , 
where the phase coherence dies, and T* w J — xt, where 
the RVB gaps appear—, and we conjecture that it is the 
crossover or "onset" identified by Ong et al. The en- 
ergy involved is relatively easy to estimate for low dop- 
ing, where gt -C <r, A. Here the average kinetic energy 
of a state with arbitrary orientation of the r's is zero, 
since there will be no correlation between kinetic energy 
and occupation; while with the locked configuration , the 
kinetic energy per electron will be of order gt. We con- 
jecture that the locking temperature T*' will be of order 
gt, so that it will rise more steeply than T c w gA . This 
agrees roughly with Ong's observations (see Fig. 2). This 
rise will stop before gt « A « ( J — gt), where the two 
contributions to £ become comparable; this will come at 
about 1/2 of the maximum doping. From this point on 
the calculation of the energy becomes quite complicated 
and we will postpone it to a later paper. Kotliar and Liu 
suggest that the onset may not change much with further 
doping 5 . 

The electrodynamics of the partially locked state is 
in principle the same as that of the vortex liquid state, 
in terms of symmetry. It has in common with that 
state that there is always an anomalous amplitude, i.e., 
a "gap" , but that the phase of this gap is fluctuating 
freely; there is a classical liquid of vortex lines. Not 
much thought has been given to the responses of such a 
liquid, particularly not to the diamagnetism, which must 
approximate that of the Abrikosov state at the melting 
line; how it diminishes from this has not been discussed. 

I have concluded that an approximation to its behavior 
may be obtained in the following way. The instantaneous 
response to an electromagnetic field is identical to that 
of an ordinary superconductor, but the current-current 
correlation decays with a finite correlation time r. From 
this picture we can predict two responses. The supercon- 
ducting response to an electric field is the acceleration 
equation: ^ = p s E ; correspondingly, for the locked but 
non-superconducting state we should have 

J s = p s rE; a s = p s r . (4) 

We have no apriori reason to select a magnitude for 1/r, 
but it certainly should remain less than A, otherwise 
the picture is meaningless; and experimental data on in- 
frared response in the RVB region suggest that the vor- 
tices freeze at 1/r = T = TjM,. 

The diamagnetic response is more interesting. Our pic- 
ture of the pseudogap- vortex liquid state is that the order 
parameter remains non-vanishing but has no long-range 
phase order. If one were to introduce a flux quantum 
instantaneously at a point one would induce supercur- 
rents throughout the sample, including the long-range 
current (h/m)p s \7(f> cx 1/r. The divergent energy due to 



such currents means that in a superconductor they must 
be screened out by a surface current, which is respon- 
sible for the magnetization of the superconductor. It is 
these magnetization currents which correspond to the ac- 
celerated supcrcurrent in an electric field and exhibit a 
delta-function response S(oj) — > l/(iu> + 1/r) when the 
gap fluctuates thermally with a correlation time r. This 
response is proportional to the magnetization which for 
low fields varies as B In ^ a . Thus the susceptibility 

X oc rln; ^ , (5) 

will have a logarithmic divergence as B — ► 0. This singu- 
larity was first observed by Ong and Lu 1 - 1 , in fact it was 
the search for a mechanism for this behavior which led to 
the present line of thinking (although the observed sin- 
gularity seems to be more precisely described by a power 
law.) They also observe that the upper critical field H C 2 
is continuous through the melting line and remains large 
and finite to temperatures where the susceptibility can 
no longer be measured. This field must be that where 
the charge becomes unlocked from the RVB, since I can't 
see how the field would have much effect on the RVB. It 
is hard to find any other explanation for the remarkable 
fact that the susceptibility remains singular over such a 
wide range of temperatures; no explanation in terms of 
critical behavior seems viable. 

It is possible to improve somewhat on the simple ex- 
pression which is meant only to show the nature of 
the singularity at B = 0. In the d-w&ve superconductor, 
there is a distribution of gaps and therefore we may ap- 
proximate its response by summing over a distribution of 
H C 2's. To get a notion of the form of the variation, let us 
assume that the distribution of the gaps is uniform; if we 
do this, we find that the expression for the magnetization 
becomes, 

M oc B (VB - 1 - In VB) , (6) 

where B = H / H C 2 max . This is illustrated in Fig. 4 and 
its resemblance to some of Ong's data is evident. 

In summary, I have proposed a theory which accounts 
for most of the anomalous experimental facts about the 
pseudogap state, and connects rather seamlessly to the 
only successful microscopic theory of the superconduct- 
ing state. Neither theory places much emphasis on the 
complications of the various inhomogeneous phases which 
tend to occur in these systems but which seem to involve 
smaller energies and weaker perturbations than the more 
striking and universal effects we discuss. (In fact, at least 
one such phenomenon, the "checkerboard" , seems to re- 
ceive a natural explanation within this theory^.) The 
crucial step seems to be the idea of visualizing the RVB 
as a separate entity in the spin Hilbert space which is 
locked to the electron charge Hilbert space by a rela- 
tively weak force. This makes irrelevant the many "ghost 
vacua" of the spin system which have led other work- 
ers astray. Yet when the spin system chooses a physi- 
cally different "vacuum" -the antiferromagnetic state, for 
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instance-the whole complex of superconductivity com- 
pletely disappears^!. On the other hand, as Ogata has 
shown, it is possible for antiferromagnetic order to coex- 
ist with d-wave superconductivity 2 ^. 

In terms of the symmetry classification of phases, on 
which often great emphasis is placed, I have not dis- 
covered any "new" entities. The superconducting state 
is only unconventional in that its quasiparticles are not 
hole-particle symmetric. The pseudogap phase is "only" 
an enormous extension of the vortex liquid. Yet each ex- 
hibits new and unexpected physical phenomena. One is 
tempted to speculate that through too much emphasis 



on mathematical categorization physical reality has been 
neglected. 
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FIG. 2: The "dome" and the region of observation of a vortex 
Nernst effect, as a function of doping for two superconducting 
systems. 
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FIG. 3: The charge and spin pseudovector triads, and the 
relative orientation when "locked" (in the superconducting 
state.) 
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FIG. 4: Hypothetical variation of diamagnetism with mag- 
netic field in pseudogap phase. Here, B — H /i/ C 2max, and the 
scale of M decreases as t(T). 



